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Abstract 

We follow up on the analysis of Mecozzi and Bellini (arXiv:1110:1253vl) where they showed, 
in principle, the possibility of superluminal propagation of neutrinos, as indicated by the recent 
OPERA result. We refine the analysis by introducing wave packets for the superposition of energy 
eigenstates and discuss the implications of their results with realistic values for the mixing and 
mass parameters in a full three neutrino mixing scenario. Our analysis shows the possibility of 
superluminal propagation of neutrino flavour in a very narrow range of neutrino parameter space. 
Simultaneously this reduces the number of observable events drastically. Therefore, the OPERA 
result cannot be explained in this frame-work. 



We dedicate this paper to the memory of Raju Raghavan who has made fundamental 
contributions in the area of neutrino physics, and who passed away while we were writing 
this paper. 
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I. INTRODUCTION 

The recent announcement of the OPERA result [l| indicating possible superluminal prop- 
agation of neutrinos has excited considerable interest. Various aspects of the experiment, 
the analysis of the data and their interpretation, must be subjected to a thorough exami- 
nation since the result has important repercussions on fundamental physics. Furthermore, 
independent confirmation or refutation by other experiments is absolutely essential. 

If the result stands, one must first see whether it can be understood within the usual 
framework of physics before giving up cherished notions such as Lorentz invariance. In 
fact, this is possible, in principle, as was shown by Mecozzi and Bellini [2|. By considering 
the interference between the different mass eigenstates of the neutrinos they showed that 
superluminal propagation is possible. We extend their analysis by explicitly including the 
effect of the finite width of wave packets and provide a numerical estimate of the effect 
for realistic neutrino parameters with three generation mixing and matter effects taken into 
account. 

Our numbers indicate that there is a very narrow region in the allowed parameter space 
with three neutrino flavour mixing in which superluminal propagation is possible in princi- 
ple. However, the survival probability for neutrinos with superluminal velocities is almost 
vanishing rendering them unobservable in practice. Furthermore, this depends crucially on 
the ratio of the distance of propagation and the energy of the neutrinos. At distances and 
energies corresponding to the OPERA experiment, superluminal propagation is not possible 
with the present limits imposed by the neutrino parameter space. If the OPERA result is 
confirmed it would require new physics. Therefore, we would like to stress the importance 
of further studies along the present lines since the OPERA experiment has opened up a 
new window on neutrino physics, which may be called neutrino optics and which should be 
pursued by future experiments. 

In Sec. II we outline the group velocity calculation in the wave packet formalism. In the 
context of neutrino oscillation, this formalism has been discussed in detail in Ref. [3|. We 
derive results for short and very long base-line propagation of neutrinos, based on which 
we present a realistic numerical analysis in the framework of 3 neutrino mixing in Sec. III. 
Other issues related to group velocity measurements will be discussed in Sec. IV while Sec. V 
concludes the paper with some remarks on the implications of OPERA-type experiments in 
future. 

II. CALCULATION OF THE GROUP VELOCITY 

The neutrino flavour states |z/q,), a = e, fi,T are related to mass eigenstates |z/j), i = 1, 2, 3 



where U is a unitary matrix. For a neutrino that starts as a flavour state « at t = 0, the 
state vector at time t is 



by 
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where we have superposed the three energy eigenstates with the same momentum p and 
then superposed different p with an amphtude of gaussian form 6"^^"^°-* ^'^ to form a wave 
packet. We have set = c = 1. The normahsation constant is determined by the condition 

{m\m) = ^- 

The probabihty amplitude for detecting lup) at time t is 



A 



Pa 



N j dpJ2Q 



ipx-iEit ^-{p-po)"^ /a? 



^ipx-iE,t g-(p-po)^/a^ 



(3) 
(4) 



where we have used 

We expand the energy E around the peak po of the gaussian and keep upto the quadratic 
terms (assuming the width of the wave packet in momentum space to be small enough): 



dE 
dp 



{P - Po) 



Id'^E 



po 



2 rfp2 



{p-poY 



(5) 



Po 



With this approximation, the p-integration in eq. (j4]) can be done to yield the result (ab- 
sorbing the constant factors into A^): 



Oi — ^Qi 



En 



d^Ei 

4 dp^ 



(6) 



po 



where XQ{t) = dEi/dp\pg t. Thus we have a superposition of 3 gaussians in a;-space with 
their centres travelling with 3 separate group velocities, 



dx^it) dEi 
dp 



dt 



(7) 



Po 



However, if the x-space gaussians are broad, that is, if a is small, then these 3 gaussians will 
interfere. To study this, let us define 







' 1 diX 




^'dx 



The integrations are straightforward and the result is 



-iAEijt\ -{v^-Vj)'^t^a'^ I 



(9) 



where 



AEjj = E^i — Eoj 



a? d^{Ei-Ej] 
4 dp'^ 



Po 
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Defining the overall-group velocity of neutrinos generated as at time t = and detected 
as 1^/3 at time t as Vj^a as was done by Mecozzi and Bellini j^, we get 



I3a 



13a ^fia* I 
3 



(10) 



This is the main result of the paper. 

We now consider the special case of two generation mixing with /i, r as the two neutrino 
flavours. Substituting a = /3 = /i, and denoting C^^ = cos^ 9, Cg'^ = sin^ 9, we have 



vi cos^ 9 + V2 sin^ 9 + {vi + V2) sin^ 9 cos^ 9 cos^AEut) g-^^'i-''^)'*''^' 
cos4 9 + sin^ 9 + 2 sin^ ^ cos2 ^ cosiAEut) e-(''i-'^2)'*''^V8 



Let us consider the two extreme limits from the above equation. First consider the case 
when the width of the wave packets in the x-space is large compared to the distance of 
separation between the two centres of the wave packets, that is a{vi — V2)t -C 1. In this limit 
eq. f ITT]) becomes 



vi cos^ 9 + V2 sin'' 6* + (f 1 + V2) sin^ 9 cos^ 9 cos(A£'i2t) 



cos^ 9 + sin^ 9 + 2 sin^ ^ cos^ 9 cos(AEi2t) 



(t^i + V2) + {vi - V2) 



cos 2^ 



1 - sin^2^ sin^(A^i2t/2) 



(12) 



This agrees with the result of Mecozzi and Bellini [2|]. 

However as the distance or time of propagation increases, the width of the wave packets 
in x-space becomes small compared to the distance of separation between the two centres of 
the wave packets, a(t>i — t>2)t ^ 1, and we get 



f 1 cos^ 9 + V2 sin^ 9 1 
cos^^ + sin^^ ~ 2 



{vi + V2) + {vi - V2)- 



cos 29 



I sin^ 29 



(13) 



which is the weighted average of the group velocities of the two wave packets. Interest- 
ingly, the effect of the gaussian suppression is precisely the same as taking the average over 
energy and distance in the denominator of eq. ( IT2l) . This makes sense since, for instance, 
at astrophysical distances, the neutrino wave-length is small compared to the distance of 
propagation. 

Thus the result in eq. f lTOj) generalises the result of Mecozzi and Bellini to wave packets of 
finite width. Because of the approximation made in eq. (|6]), for large a the damping factor 
exp(— (fi — V2)'^t^a? /9>) is only approximate, although its exact replacement also will damp 
the oscillatory factor cos(Ai?i2t). 

If the calculations are done for waves of infinite spatial extent, the integrals occurring 
in the numerator and denominator of eq. (|H]) would be individually divergent, although the 
final result would turn out to be finite. However, it is better to do these calculations with 
wave packets of finite width a, as has been done above, and take the limit of a —t- in the 
end. 

We have superposed the 3 mass eigenstates with the same momentum p but different Ei 
to form the neutrino state vector in eq. (j2]). Should one superpose different pi but the same 
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E, or, different pi and different Ei7 This question has been studied in recent hterature 
Such possibihties will be included in a future study of the group velocities of neutrino waves, 
which is under preparation. 

For neutrinos whose energy is very large compared to their rest masses, the formula in 
eq. (fTTi) may be written in the form 



f 1 cos"^ 9 + V2 sin^ 9 + {vi + V2) sin^ 9 cos^ 9 cos y e 



cos^ 9 + sin^ 9 + 2 sin^ 9 cos^ 9 cos y e-^^vlv)WI^ ' ^ ' 

where y = A'm'^L/(2E), Am^ and L being the measured mass-squared difference and base- 
line distance and 6p = a/y/2 the width of the p-space wave packet. Thus the exponential 
damping factor multiplying the oscillatory factor cosy is simply e~^^^^^^ ^ where Sp/p is 
the fractional uncertainty in momentum. 

We also note that the formula in eq. fllip and eq. (Iiup for two- and three- generations are 
valid, to the leading order, for propagation through matter at constant density when vacuum 
values of Ei and the mixing angles are replaced by their matter-dependent values. 

From eq. ([6]), the normalised oscillation probability, that is, the probability for detecting 
flavour /3 at time t is 



I3a 



I \Ap^\^dx = J2 |Cf"r + $^2Re (cf"cf *e-^^^-*)e-(^^/^')'3''/^ . (15) 

i i>j 



This differs from the usual oscillation formula by the factor exp{ — {6p/ p)'^y'^ / 4) in the second 
term. In view of the successful neutrino oscillation phenomenology achived so far, we will 
assume that a = \/2Sp is so small that this factor can be replaced by unity for all the 
terrestrial experiments as well as solar neutrino experiments. 

Since the oscillation probability given in eq. ( ITSl) is the denominator in eq. (fTOll . the 
group velocity can become very large if the oscillation probability is very small. In fact, it 
can become infinite if the oscillation probability at that distance is zero. This is the origin 
of superluminal propagation, as our analysis in the next section will clearly show. 

We now come back to the interpretation of {x)i3a defined in eq.(IH]) which is the basis of the 
above formulae. Note that in the denominator we have / (ixlA^Qp instead of Yl/s J dx\Ai3a\'^ 
which is unity. Thus (x)/3q must be interpreted as the conditional measurement of the position 
of the neutrino under the condition that only Ua is detected. Here the probability amplitude 
for detecting it as Ua itself is regarded as the wave function for normalising the expectation 
value of X. This is to be contrasted with the usual definition of the expectation value of x, 
independent of the flavour detected, namely 



/3 



We may distinguish this case by calling it the unconditional measurement of the expectation 
value of X. 

Before we go to the numerical analysis, we make some general remarks. As already 
pointed out, the origin of the superluminal propagation is the vanishing of the oscillation 
probability Pjsa in the denominator of eq. (ITOl) . In other words, superluminal propagation 
and the vanishing of the oscillation probability go together. Since the number of events also 
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vanishes, one has the paradoxical situation of unobservable superluminal propagation. Our 
numerical analysis below is subject to this criticism. 

In the particular case of the OPERA experiment, P^a becomes the survival probability 
Pnn. There is no evidence for the survival probability P^^ in OPERA becoming vanishingly 
small. Hence explanation of superluminality through enhancement of the group velocity in 
the standard oscillation frame-work as discussed above is untenable. 

Our numerical analysis will be based on the group velocity derived from the conditional 
measurement of (x)/3q, defined in eq. ([H]), following Mecozzi and Bellini j2|. Alternatively, one 
could base the analysis on the group velocity derived from the unconditional measurement 
{x)a defined in eq. (ITB]) . Since this does not have the vanishing denominator, it will not lead 
to superluminal group velocity. 

Actually, for the OPERA experiement, both (x)/?^ and {x)a will give essentially the same 
result for the group velocity since P^^ at OPERA does not deviate very much from unity, 
as far as is known. 



III. NUMERICAL ANALYSIS WITH THREE GENERATIONS 



To obtain a realistic estimate for the group velocity of muon neutrinos, and its implica- 
tions, we discuss the three generation scenario. 

As shown in the previous section, the gaussian smearing will have no effect if {Sp/pY <^ 1. 
We will work under this assumption and make estimates for this quantity later, showing that 
it is indeed small. As a first approximation we neglect the matter effects. 

Following eq. flTU]) . we may write the "group velocity" of the superposition that starts as 

and is detected as (which we denote as for simplicity) as 



= V2 + Av 



The factors Sij are given as 
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where the denominator, which is simply the survival probability of i/^, is given by, 

3 

j>i=i 

Here 

AE.A (E. - EM. 

1.27- 



AE,,t _ {E, - Ej)t _ ^ ^ A,,(eV2)L(km) 



E(GeV) 



(17) 

(18) 
(19) 

(20) 
(21) 



where Aij = mf 

lepton mass matrix is diagonal as 

C13C12 



m| and the mixing parameters are given in the basis where the charged 



U 



C13S12 Si3 e 

-C23S12 - S23S13C12 e'^ C23C12 - S23S13S12 e'^ S23C13 
■S23'5i2 ~ C23S13C12 e* — S23C12 — C23S13S12 e* C23C13 



i8 



(22) 
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FIG. 1: (L) The terms and 832 as a function of 623- The values on the left are for the normal 
hierarchy solution with ^13 = while those on the right are for the inverted hierarchy with ^13 = 10°. 
In both sets, the curves for 5*32 have been offset by 623 = 0.1° for clarity, else the two curves for 
S12 and 5*32 would overlap each other. In the former case, 5cp = and L/E = 611.165 km/GeV 
while for the latter, dcp = 180° and L/E = 604.629 km/GeV. (R) The denominator P^^ plotted 
as a function of 623 for the same parameters. The dips as — )■ correspond to the sharp peaks 
in Sij on the left. 

where Cij = cos 6ij, Sij = cos 6ij and 6 is the CP phase. The coefficients f/^j correspond to 
the second row of the mixing matrix. Note that in the absence of any mixing in eq. (fTTj) 
reduces to V2, as it should. 

We now present a numerical analysis of the possible values of Vf^ with realistic parameter 
values for the mixing angles and mass squared differences. We use typical values/ranges 
for the mass-squared differences and neutrino mixing parameters: A21 = 7.6 x 10~^ eV^, 
IA32I = 2 X 10^3 eV^, 012 = 34°, 36° < 623 < 54°, ^13 < 10° while Scp is unknown Q. The 
sign of A32 is not known and there are two possible hierarchies, m| > known as normal 
hierarchy (N) and m| < known as inverted hierarchy (I). We consider results for both 
these cases. 

Note that the propagation can become superluminal when the Sij are significantly en- 
hanced over unity. Fig. [1] shows the values of 6'23, close to maximal mixing, where this 
enhancement is seen, for different mass hierarchies as well as values of ^13. (Notice therefore 
that this measurement is sensitive to the mass ordering in the 2-3 sector which is as yet 
unknown) . 

This enhancement of the group velocity occurs for L/E (km/GeV) ~ 600, well within the 
range of the OPERA experiment, viz., 7 < L/E{ km/GeV) < 730. 

However, it can be seen from the right hand panel of Fig. [1] that P^^ is very close to zero 
precisely at these values, that is, the enhancement of Sij is dominantly due to the vanishing 
of the denominator. (The maxima of Sij are slightly offset from the minima of P^^ because 
of the dependence of the numerator on the mixing parameters as well. In fact, when ^13 = 0, 
Pfj_^ and Sij both vanish near 623 = 45°.) This in turn means that there will be hardly any 
events that survive at these values, for this conditional measurement! 

Several remarks are in order: 

• The first is that the mixing parameters as well a.s L/E need to be extremely fine-tuned 
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in order to get this effect. However they are almost washed out by the suppression of 
the event rate due to small survival probability. 

• The extreme enhancement of Sij over unity is required in order to obtain superluminal 
velocities commensurate with the OPERA observation. This is because the coefficients 
of Sij which are (vi — Vj) in eq. ( IT7I) behave as A^j/ (2p^) and are small due to the highly 
suppressing factor of in the denominator. It is unlikely that such small excesses 
will be measurable by any experiment. The larger the enhancement required, the more 
finely tuned are the parameters. 

• Moreover, the fine-tuning in the value oi L/E would imply that the enhancement of 
the velocity only occurs for a single E value when the base-line distance L is kept fixed, 
in contrast to the observed roughly constant enhancement over a range of 1 < £^ < 100 
GeV as observed in OPERA 

One way of working around this limitation is to consider the energy dependence of the 
matter-dependent contributions. The usual electro-weak interactions in matter lead to 
the inclusion of a matter-dependent potential, Vew, that alters the matter dependent 
mass-squared differences in a non-trivial way; however, the resulting change is not large 
enough to give the required enhancement of the factor of 14 in the term A^L/(4i?) at 
{E) = 17 GeV. 

One possible solution is the inclusion of a new matter interaction Vnew with the same 
energy behaviour, but about an order of magnitude larger than the usual electro-weak 
potential. Then, at energies around (E) = 17 GeV relevant to the OPERA data, the 
expression for the matter-dependent mass-squared differences simplifies to 

A™ ^ Aij + 2EVnew , (23) 

so that the term occurring in Av^ becomes 

AE ^ AE 2 ' ^ ' 

The energy dependence of this matter potential is exactly what is required so that 
the ratio A'^/E is approximately independent of E for E ~ 10s of GeV (but does 
not significantly alter the atmospheric neutrino analysis). Such an energy-independent 
contribution would remove the fine-tuning in L/E that currently occurs in the expres- 
sions for Av^ and would in principle uniformly allow for superluminal propagation of 
all velocities relevant to the OPERA analysis. However, while resolving the fine-tuning 
in L/E, the new matter potential enhances the relevant terms in a manner identical to 
the vacuum analysis, viz., through the vanishing of P^^. Hence, it runs into the same 
difficulties with observing the effect as discussed earlier. 

• In spite of what is said above, superluminal group velocity is a real effect and may 
be observable in future, as shown in the following example. Consider the 2-flavour 
case where we neglect m2 and set = A32 ~ 2 x 10~^ eV^. For smaller energies of 
the order of KeVs, observable excesses of over unity can be obtained for a modest 
value of P^^ ~ 0.03. This corresponds to sin^ 2^23 = 0.97, with the ratio L/E tuned 
so that sin^(A32L/(4£')) = 1 (achievable with energies in KeV and length in meters); 
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this leads to t^^ = 1 - A32/(2p2) x (0.5 - cos 2^23/^^^^) ~ 1 + 6 A32/(2p2). For energy E 
in KeV we get an enhancement v^ — l = Q 10~^/-E^; it may thus be possible to observe 
the effect since P^^ is not zero. Of course, there are many experimental problems to 
be overcome before such an observation is achieved. 



IV. OTHER IMPLICATIONS OF GROUP VELOCITY MEASUREMENT 

We add a few remarks on the implication of the OPERA result independent of the fact 
that it is superluminal or subluminal. Consider a possible measurement of group velocities 
of electron and muon neutrinos in a future possible experiment. Following eq. ( IT7|) . these 
velocities in vacuum may be written in the form 

Ve = Vi + S2i{v2 - Vi) + 5*31 (t;3 - Vi) ] 

Vf, = V2 + 5'f2K - V2) + 5*^2(^3 - V2) , (25) 
where S^j are given in eqs. (ITSl) and (HM and S^j are given by 

^21 = : ,rr 101.. 10 • 2.. ^ , /^X ' (26) 

and 
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\^sin\AEijt/2) 



^ l^esHl - 2|^eipsin^(A^i3t/2) - 2 1 [7^2 P siu^ ( AE23t/2)] 
' 4E->,=i|f/e.P|t/e,f sin2(AE,,t/2) 



Simultaneous measurement of these two velocities immediately gives information on the 
ordering of masses m2 and since we already know that mo > mi from the solution to the 
solar neutrino problem. Unlike earlier proposed solutions |8|, this does not require matter 
effects to resolve the issue. Even a short base-line experiment with neutrino factories with 
muon storage rings may help resolve the hierarchy issue. Note that the Ve and could refer 
to either neutrinos or antineutrinos. The caveat is that the effect is magnified only for a 
given set of parameter values including L/E. 



V. CONCLUSIONS AND SOME REMARKS 

To summarise, we have considered the superposition of 3 energy (mass) eigenstates with 
the same momentum p, but with different energies Ei to form a neutrino flavour state. We 
have included the effect of the finite width of wave packets of the same width and considered 
its effect on the group velocity of the neutrino flavour state. 

When this is applied to the propagation of for distances of the order of hundreds of kms 
and energies corresponding to the observation of muon neutrinos in the OPERA experiment 
we find the effect is very small. Effectively the finite width can be ignored as is done in 
Ref. 0. 

However the effect of width increases with distance. For astronomical distances the effect 
of finite width is to reduce the group velocity to a weighted average of the velocities of 
the individual energy eigenstates. This explains why, even if the OPERA result is correct. 
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there is no contradiction with the absence of superluminal propagation of neutrinos from 
supernova SN1987a. 

Our analysis shows that superluminal propagation of neutrinos occurs whenever the oscil- 
lation probability corresponding to that measurement vanishes. For instance, muon neutri- 
nos are dominantly observed at OPERA; these will exhibit superluminal behaviour precisely 
when the survival probability P^^ vanishes. While P^^ does not vanish for the parameter 
range of the OPERA measurement, in principle, this feature makes superluminal neutrinos 
effectively unobservable through conditional measurement as explained in section 2. 

Finally we discuss the implication of the OPERA result irrespective of superluminality 
or otherwise. The neutrino mass ordering or hierarchy is not fully known. The present 
understanding is that this requires separate measurement of oscillations of neutrinos and 
anti-neutrinos in the presence of matter. The OPERA measurement of neutrino flavour 
"velocity" has added a new way of precision measurement of neutrino parameters and may 
have significance in the context of neutrino mass hierarchy (provided the parameters lie in 
an extremely narrow part of the known range of neutrino parameters). A detailed analysis 
including the effect of matter will be presented elsewhere. 

After most of the work reported here was completed, we came across Ref. [9| which overlaps 
with parts of the present paper. We thank Tim R. Morris for bringing to our attention his 



paper Ref. [10| in which similar ideas have been discussed. In particular, the extreme fine- 
tuning of parameters required to produce velocities of the order observed in OPERA and 
the difficulty that the effect is seen only at fixed energy was pointed out in this paper. This 
paper points out the existence of multiple peaks, that would give rise to a strong energy 
dependence, in addition. 

We are grateful to G. Date for many discussions, clarifications and comments. We also 
thank N.D. Hari Dass, S. Kalyana Rama, D. Sahoo, N. Sinha and R. Sinha for many discus- 
sions. 
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